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The problem of non-perturbative description of unsteady premixed flames with arbitrary gas 
expansion is solved in the two-dimensional case. Considering the flame as a surface of discontinuity 
. . . with arbitrary local burning rate and gas velocity jumps given on it, we show that the front dynamics 

OO ' can be determined without having to solve the flow equations in the bulk. On the basis of the 

, Thomson circulation theorem, an implicit integral representation of the gas velocity downstream 

■ is constructed. It is then simplified by a successive stripping of the potential contributions to 

04 ' obtain an explicit expression for the vortex component near the flame front. We prove that the 

(— I ^ unknown potential component is left bounded and divergence- free by this procedure, and hence 

. can be eliminated using the dispersion relation for its on-shell value (i.e., the value along the flame 

front). The resulting system of integro-differential equations relates the on-shell fuel velocity and the 
\ front position. As limiting cases, these equations contain all theoretical results on flame dynamics 

, established so far, including the linear equation describing the Darrieus-Landau instability of planar 

flames, and the nonlinear Sivashinsky-Clavin equation for flames with weak gas expansion. 

Ch ! PACS numbers: 47.20.-k, 47.32.-y, 82.33.Vx 
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I. INTRODUCTION 



Although the elementary physical mechanisms underlying flame propagation are well-understood by now, a global 
mathematical description of the process is extremely difficult. The reasons - which flames share with condensation 
discontinuities in supersaturated vapors [l| , ablation fronts driven by lasers or X-rays , thermonuclear waves in type- 
^ i' la supernovas Q or rapid decomposition of explosive liquids [l| - can be summarized briefly: all involve propagating 
free-boundaries crossed by nonzero mass fluxes and separating subsonic flows of markedly different densities. In 
_ -, ■ view of this it is not surprising that, since their identification as self-propagating slow deflagrations (as opposed to 

CN ■ „„N .1 _ ioor,„ nil a U.J ac^ l ^ • 



detonations) in the 1880s [1], flames had to wait some flfty years before their theoretical study began at the simplest 
level of linear stability analysis of planar flames 0, Q . It took some forty years more for the first consistent account of 
p\| , nonlinear effects to appear. It was shown by 0] how these effects can be described in the case when the fresh to burnt 
' gas density ratio, 6, is close to unity. The latter condition is a principle limitation for a perturbative treatment of 
[ nonlinear saturation phenomena. Indeed, the weak nonlinearity expansions arc simply self-contradictory in the case 
[ — ■ of steady flames with a finite {d —1) [1]. Despite the fact that the practically relevant values of are 5 to 8, the small 

' {6 — 1) approximation has so far been the sole theoretical method^ available to handle the flame front dynamics. 
C , One of the essential difficulties encountered in any analytic treatment of flames is the virtual impossibility to solve 
■ the flow equations governing the dynamics of the burned gases. This is an exceedingly complicated problem which 
^ . requires flnding solutions to a system of nonlinear partial differential equations in the regions ahead of the flame front 
and behind it, to be chosen so as to satisfy a number of jump conditions expressing the conservation of mass, energy 
. , and momentum across the moving front. The front dynamics itself is determined by the so-called evolution equation 
' descri bing the local fuel consumption rate as a functional of the fuel velocity distribution along the front and the front 
■ ■ ' shape [lO|. Even if the gas flow is potential upstream, as is the case for flames propagating in an initially quiescent 
fluid, this property is lost in the downstream region because vorticity is generated by the curved flame front, so that 
the problem of solving the flow equations is faced in its full generality. 

Struggling with this problem is indeed unavoidable if one is interested in the explicit structure of the burnt gas 
flow. However, it is the evolution of the flame front, its position and shape that usually constitute the main concern 
in practice. This limitation of the problem raises naturally the following dilemma [see |Tl|, [l^: On the one hand. 



We do not consider models inconsistent from the standpoint of the fundamental equations, such as for instance the Frankel's Q potential- 
flow model. The latter assumes finite values of {6 — 1) and zero vorticity production in the flame front at the same time; these two 
assumptions arc only consistent in the lowest order of the small (8 — 1) expansion, in which case the Frankel equation reduces to the 
Sivashinsky equation 0. 
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deflagration is an essentially non-local and nonlinear process with all its complications mentioned above; on the 
other hand, this non-locality itself is determined by the flame front configuration and the gas velocity distribution 
along it, which play the role of boundary conditions for the flow equations and thus control the bulk flow. In such 
circumstances, is it really necessary to know explicitly the flow structure in the entire downstream region in order to 
describe the front evolution in a closed form, i.e. in a form involving only quantities defined on the flame front? 

As to the steady configurations, this question was settled out in the negative by [111, More precisely, it was 
shown that the only piece of information about the gas flow downstream that is really necessary to derive an equation 
for the flame front position is the value of the vortex component of burnt gas along the front or, using the terminology 
used by [l^ . the on-shell value of this component. The remaining unknown potential component of the gas velocity 
is excluded eventually using a "dispersion relation" for its on-shell values (expressing analyticity of this component 
in the downstream region), thus providing one with a description of flame shapes in a form which is closed in the 
above-mentioned sense. The purpose of the present paper is to generalize this construction to the case of unsteady 
flame propagation. Although the procedure is essentially the same as in the steady case, a subtle point is worth to be 
emphasized. According to [llill^l, the derivative of the vortex component along the flame front is local, i.e., its value 
at a given point is a function of the on-shell fuel velocity, its derivatives, and the front shape at the same point. In 
view of this, one might expect the generalization to the unsteady case to be purely kinematic, namely, that it would 
amount to rewriting the steady equation in terms of the gas velocities relative to the local front velocity. Wc will see, 
however, that this is not so because of a peculiar - yet unavoidable - spatial non-locality of the vortex component, 
which appears naturally in the unsteady case to account for the effect of the flame history on its current evolution. 
Because of this complication, one has to be more careful with the spatial integrations involved in the definition of 
the vorticity mode. To define the improper integrals along the flame front, we use an intermediate regularization. 
Specifically, we introduce an exponential damping of the contributions coming from remote parts of the front. This 
regularization is eventually removed using analytic continuation to the case of zero damping. 

The paper is organized as follows. We first construct an implicit integral of the flow equations downstream on the 
basis of the Thomson circulation theorem, which expresses the gas velocity in terms of its boundary values and vorticity 
distribution behind the front. This is done in Sec. [Til The integral representation is then used in Sec. IIIII to obtain 
an expression for the vortex mode of the gas velocity near the front, which is accomplished by successive stripping 
burnt gas velocity of potential contributions. We prove that the unknown potential component is left bounded and 
divergence-free by this procedure. Hence, it can be eliminated using the dispersion relation for its on-shell value, thus 
leading to the main integro-diffcrcntial equation written down in Sec. IIVI This equation relates the on-shcU value of 
the fuel velocity and the flame front position, and together with the evolution equation constitutes the closed system 
for these quantities. Finally, it is verified in Sec. |V] that the derived equation contains as simple limiting cases all 
known theoretical results on fiame dynamics, namely, the linear equation describing Darrieus-Landau instability of 
planar flames 0] , including the case when the flame propagates in a time-dependent gravitational field [Tsj , the 
nonlinear Sivashinsky-Clavin equation for flames with weak gas expansion 0, fl^ , and the stationary equation derived 
by [ll]- Appendix contains a consistency check for the results obtained in Sec. IIIII 

II. INTEGRAL REPRESENTATION OF THE FLOW EQUATIONS 

Consider a flame propagating in an initially uniform premixed ideal fluid. Our analysis below relies substantially 
on the well-known Thomson theorem stating that circulation of the gas velocity over any closed material contour 
drawn in an ideal fluid is conserved as it is convectcd. This statement takes an especially simple form in the case of 
two-dimensional (2D) incompressible flows, since not only the circulation itself, but also the value of vorticity carried 
by any fluid element is then conserved. Since the space dimensionality is not that important in the formulation (not 
the resolution) of the dilemma mentioned in the Introduction, we will be concerned in what follows with the simpler 
2D case. We will further specify our analysis to flames propagating in a channel of constant width b. Let the Cartesian 
coordinates {x, y) be chosen so that the y-axis is parallel to the tube walls, y = — oo being in the fresh fuel. These 
coordinates will be measured in units of the channel width, while the fluid velocity, v = (wi,f2), in units of the 
velocity of a plane flame front relative to the fuel. It will be sometimes useful to denote the Cartesian components of 
V by (w,u). Finally, the fluid density will be normalized on the fuel density, 9 > 1 denoting its ratio to that of the 
burnt gas. 

It will be more appropriate for our purposes to reformulate the problem under consideration as a problem of 
propagation of an unbounded spatially-periodic flame. Namely, assuming the channel walls ideal, given a flame 
conflguration described by the functions f{x,t),v(x,y,t), x G [0,-1-1], where f{x,t) denotes the flame front position 
at time instant i, using the boundary conditions /' = 0, w = for x = 0,1, we continue this configuration to the 
domain x £ [—1,0] according to 



f{x,t)^f{-x,t), w{x,y,t) = -w{~x,y,t) , u{x,y,t)=u{~x,y,t), 



(1) 
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and then periodically continue it to the whole x-axis. Note that having imposed the boundary condition /'(O) = 
/'(I) = 0, we thereby exclude the possibility of stagnation zone formation near the end points of the flame front 
(see [m for detail). We also assume that the flame is stable with respect to the short wavelength perturbations i.e., 
that there is a short wavelength cutoff, Ac. This cutoff ensures smoothness of the functions under consideration. In 
particular, it prevents the development of singularities of the front shape such as the edge points which would occur 
otherwise leading to discontinuities in the values of the flow variables or their derivatives. That Ac often exceeds 
the actual thickness of the flame preheat zone significantly [T6j has yet another virtue: the Reynolds number based 
on Ac and the fuel properties is typically over ~ 10^, and hence is fairly large when based upon the width (> Ac or 
3> Ac) of the channel where the flame studied below is meant to propagate. It then makes sense to model the flame 
as a surface (or line in 2-D) equipped with a local Ac-dependent propagation law, and embedded in ideal fluid flows. 
We shall return to this issue in the final section of the paper, merely mentioning here that viscosity effects are known 
from direct numerical simulations to have negligible influence on the shape and the speed of steady curved flames. 
In our formulation, the flow velocity obeys the following equations in the bulk 



dx. 



0, (2) 



d(J dcr 

+ (3) 

where {xi,X2) = {x,y), 

is the vorticity, and summation over repeated indices is understood. The first equation expresses continuity of the 
incompressible flows, while the other the Thomson theorem. It is not difficult to see that v satisfying Eqs. ([2]), (|4]) 
can be written in the following integral form (Cf. derivation of Eq. (9) in Kazakov[T3) 

Vi = Sikdk J dli simVm — ~ di j dlk Vk— ~ Bikdk j ds , (5) 

A A E 

where Eik = —£ki, £12 = +1, dt = d/dxt. E and A denote any part of the downstream region and its boundary, 
respectively, r is the distance between an infinitesimal fluid element ds at the point (i, y) and the point of observation 
X = (x, y) G E, = {xi — x^)^, and dli is the line element normal to A and directed outwards of E. Indeed, to evaluate 
divergence of the right hand side of Eq. ([5]), we use the relations diSikdk = 0, 

dllnr ^ 2ttS{x — x)6{y — y) , (6) 

where 5{x) is the Dirac delta-function. It follows that diVi = for any point x inside S. Similarly, evaluating the curl 
of Eq. ([5]) with the help of Eq. ^ and eik£im = 5km gives the identity SikdiVk = <J. We are now to employ Eq. ([3]) to 
rewrite the last term in Eq. Q as an integral over fluid particle trajectories. 

To this end, we specify that S is spanned by fluid elements that crossed the flame front between two points on it 
with fixed abscissas ii = — A and X2 = -'rA during the time interval [—T, t], where t is the given time instant the flow 
is observed at (see Fig.[T]). Here A^T are some positive parameters tending eventually to infinity, so that S will then 
fill the whole downstream region: E = {.t, y : y > /(i, t)}. But for the time being, we keep them finite. The improper 
integral over an infinite E will be defined later, in Sec. IIII Al Then we have, by virtue of the vorticity conservation. 



Jds a In 1 

E 



+A t 

^di j dTN[i,T)vl{i,T)a+{S:,T)\n{[x ~ X{i,t,T)f + [y ^Y{i,t,T)fY''^ , (7) 

-A -T 



where iV = I -t- {df /dxf' , v"^ = Vi+rn is the normal burnt gas velocity relative to the flame front, 

( - ^ - I ^\ / 4\ 9f{x,t) 
v+^(w+,u+), u+{x,t) = u+(x,t) — — , 
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Hi is the unit vector normal to the front (pointing to the burnt matter), and {X{x,t,T), Y{x,t,T)) is the current 
position of a fluid element that crossed the point (i, /(i, r)) on the flame front at r. It is taken into account in Eq. ([7]) 
that the "volume" ds of this element is conserved in view of the flow incompressibility, and hence can be written 
as dxdrN^x, t)v^{x, t). Changing the integration variable, t ^ t — r, in the expression ([7]), and substituting it into 
Eq. (O gives 

1 1 

Vi = Eikdk j dh eimVm-^ - J dlkVk^-^dk J dxK{x,y,x,t) , (8) 
A A -A 

where the integral kernel K is defined by 

T+t 

K{x,y,x,t) = — dTM{x,t — t) 
J 


X In {[x - X{i,t,t - T)f + [y ~Y{x,t,t ~ T)]^y^^ , (9) 
M(i,r) = N{i,t)vl{i,t)a+{x,t) . (10) 

III. NEAR-THE-FRONT STRUCTURE OF THE VORTEX MODE 

The integral representation ([5]) of the flow velocity downstream is to be used below to obtain an expression for its 
vortex component near the flame front. More precisely, we are going to define a vortex component in a way that 
would allow an explicit expression for its on-shell value in terms of the on-shell gas velocity Vi+. For this purpose, we 
decompose the velocity field as 

+ i = l,2, 

where is a potential component satisfying the following requirements (DvP denotes any first-order spatial derivative 

of vpy. 

a) div DvP = 0, 

b) rotDvP = 0, 

c) DvP is bounded, in the sense that it remains finite in the limit A,T ^ +oo, i.e., for an infinitely expanding S 
region. 

The on-shell expression for the vorticity contribution will be obtained by a step-by-stcp simplification of Eq. Q 
throwing away potential fields fulfilling the above conditions a)-c). Although this derivation follows closely that of 
p^ . we give it here in detail to make clear the point where non-stationarity of the problem comes into play. 

Let the equality of two functions ipi{x, y), ip2{x, y) up to a fleld satisfying a) - c) be denoted by ipi = Lpi. First of 
all, as we saw in the preceding section, the flrst two terms on the right hand side of Eq. ([5]) have vanishing curl and 
divergence, and hence satisfy also the requirements a) and b). Furthermore, a simple power counting shows that c) 
is also met. Indeed, consider part A\F of the contour A, where F denotes the flame front. Representing this part 
as a semicircle with radius R oo, we note that D^lnr = 0{\/E?) for any given a; e S {D^ denotes any second 
spatial derivative). Taking into account also that v = 0(1), dl = Rdcf), where 4> G (0, tt) is the angular coordinate of 
the point x on the semicircle, one sees that after spatial differentiation, the two integrals over A\F on the right of 
Eq. ([5]) vanish in the limit i? — > oo. Similar consideration shows that the same integrals over F are convergent, thus 
proving that DvP remains bounded downstream in the limit A, T — > oo. So we can write 

+A 

Vi = — ^^fe / dxK{x,y,x,t) . (11) 



Next, we note that since we are eventually interested in the on-shell value of the vortex component, we may take the 
observation point (a;,y) as close to the flame front as we like, i.e., y ~ f{x,t), [y > f{x,t)]. The vortex component at 
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such points is determined by a contribution coming from the integration over small r and x ss x. Indeed, taking the 
curl of Eq. pT|) . and using Eq. ^ yields 

+A T+t 

SikdiVk = J dx J dTM{i,t — t)S{x — X[x,t,t — T])S{y — Y[x,t,t — t]) , 

-A 

which explicitly shows that a nonzero contribution to the vorticity comes only from the point {x, r} obeying the 
equations 

X- X[i,t,t-T]=0, y -Y[i,t,t-T]=0, (12) 
which for y w f{x, t) state that the point {X, Y) is close to the flame front, and hence 

X{x,t,t-T) = X + W+{x,t)T + 0{t'^), 

Y{S:,t,t-T) = f{S;,t-T)+u+{x,t)T + 0{T^). 

Expanding also f{x,t — r) to the first order in r and omitting the symbols O(t^), we thus have the following 
approximate expression for the fluid particle trajectory 

X{x,t,t — t) ^ X + ■w+{x,t)T , Y{x,t,t ~ t) ^ f{x,t) + u+{x,t)T . (13) 

It follows that if these expressions are used instead of the exact ones to calculate the kernel K, we still have the 
true distribution of vorticity near the flame front. Indeed, it was just shown that any integration over values of 
{i,T} not satisfying Eq. where Eq. is not valid either, gives rise to a potential contribution. In particular, 
the property b) of this contribution is preserved, and so is the property a), since the right hand side of Eq. (fTTj) is 
divergence-free identically whatever the form of the kernel. Finally, it is not difficult to see that the condition c) is also 
satisfied. Indeed, the above transformation of exact trajectories into the straightened ones given by Eq. p3p leaves 
the expression (|lip bounded. Therefore, the potential field added in the course of this transformation is bounded, 
too. It should be stressed that from now on we are concerned only with the on-shell value of the vortex component, 
so the meaning of the symbol = should be further specified. Of course, the transformation of trajectories changes 
the bulk vorticity distribution, and thereby the velocity field downstream. However, after the transformation as well 
as before it, integration over finite rs in Eq. ([9]) results in a field that is potential near the flame front. Both fields 
are of complicated structure which is unknown in general, but since they are potential near the front and bounded, 
we can use the on-shell value of their difference to define a field satisfying a)-c) in the whole downstream region. 
The existence of this field is guaranteed by the Cauchy theorem. Namely, we use the Cauchy formula to construct 
the field satisfying a)-c) as the analytic function with the given boundary value. It is in this sense that the above 
transformation is said to respect the properties a)-c). In particular, the symbol = is used below to relate the on-shell 
values (or near-the-front values, in the case of the integral kernel) of functions that differ by a field satisfying a)-c) 
downstream. 

We now proceed to an explicit evaluation of the integral kernel ^ which, after the transformations (|13p are 
performed, takes the form^ 

T+t 

K{x,y,£,t) = - J dTM{S:,t-T)\n{vlT^ -2{r ■v+)T + r^Y^'^ , (14) 



where r = (x — x,y — /(i, t)). The integrand here involves M(i, t — t) which is an unknown function of r. In view of 
what has been said about the near-the-front structure of the vorticity mode, one might think that it would be sufficient 
to expand this function to the first order in r, and then calculate the integral. However, this operation is not allowed 
as it would violate the condition c) and, as a result, would yield erroneous predictions (see Sec. 5). In particular, the 
on-shell value of the vorticity component cannot be found by setting r = in the argument of M. There is no such 
problem in the case of steadily propagating flames, as M is then time- independent in a frame of reference attached 



^ When transforming the kernel K{x,y,x,t), we use the same symbol = to relate expressions that upon substitution in Eq. (lilt give rise 
to fields that arc equal in the sense of = . 
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to the front. To overcome this difficulty, we will explicitly extract the singular part of the t- integral, which is related 
to the singularities of the logarithm located at the points 

T± = — (n±iVl-^A , n = ^^^—^ (15) 

in the complex r-plane. For this purpose, we first formally integrate Eq. (|14p by parts 

T+t I T+t \ 

K{x,y,x,t) = j di j dTiM{S:,t-Ti)\\ii{vlT^ -2{r-v+)T + r^Y^^ 

T+t 

Inr f , , 
= / dTiM{x,t- Tl) 



T+t / T+t 

^ ^ dr i [ dTiM{x,t- Ti) 



2tT J \ J ;[r — T_|_ T — T 

\ T J 

The first term on the right gives rise to a pure potential which is bounded, so it can be omitted. Denoting also 

T+t 

rfTiM(i,i-ri) = M(J,r,i), (16) 



/ 



we thus have 

T+t 

K{x,y,x,t)^^ J dTM{x,T,t) 



To extract the singular part of this integral, we deform the contour of integration in the complex r-plane so as to 
move it away from the poles; this of course implies certain time- wise restrictions on the function M(i, r, t), and hence 
on M{x,t). We shall return to this matter later (see Sec. 7 and Appendix). Here we mention only that in essence, the 
function M (i, t) is required to be analytic in a vicinity of the real axis in the complex i-plane, which is guaranteed by 
the existence of nonzero A^. To respect reality of the kernel, we take the singular part as a half sum of two expressions 
obtained respectively by deforming the contour above and below the real axis. The contribution we are interested in 
comes from moving the contours beyond the poles [see Fig. [2] . 
According to the Cauchy theorem, 

K{x,y,x,t) = {M(i,r+,i) - M(x,r_,t)} 

47r 



.1. [ dTM{x,T,t)l^— + ^ 



C-UC+ 



(17) 



Instead of proving that the integral over the contour C_ U C+ is free of singularity and gives rise to a bounded 
divergence-free potential, it is easier to show that the first term in Eq. (|17p is bounded and reproduces correctly the 
vorticity distribution at the flame front. This is done in Appendix. Thus, 

K{x,y,x,t) = I dTM{i,t-T). (18) 
Inserting Eq. p8)) into Eq. (fTTj) yields finally 

+A -r+ 

< = l£^k^k J dx J drMii, t-T). (19) 
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Since rjji = r_ , the right hand side of this relation is real. 

To make the meaning of the calculation performed more vivid qualitatively, it is useful to mention an interrelation 
between the roles the conditions a)-c) play in the above derivation. When calculating the near-the-front value of the 
kernel K(x,y,x,t), we retain terms of the first order with respect to r. This is sufficient for the calculation of the 
vortex component of velocity at the front, taking into account that this quantity is determined by the first spatial 
derivatives of the kernel, and that r = tv+ at the point defined by Eq. ()12p . This means, in particular, that t 
in the function M{x,t — r) cannot be neglected. As was mentioned above, this function cannot be expanded in r 
either without violating c): the condition M{x,t) — > for t — > — oo guarantees convergence of the r-intcgral in the 
limit T oo. At the same time, it is seen from Eq. (fT9|) that as a result of the r-integration, dependence of the 
function M{x,t — r) on r is transmuted into coordinate dependence. This dependence does not affect the vorticity 
distribution along the front, because r = r/v^ along the streamlines, so that t = when the observation point is 
taken on the front (r = 0). Thus, the seemingly innocent condition c) entails a nontrivial change in the structure 
of the potential component of the burnt gas velocity in comparison with the steady regime. In the latter case, on 
the other hand, the condition M{x,t) ioi t ^ — oo docs not apply, but since the M-function is independent of 
time, it remains independent of the coordinates [x, y) at all stages of the calculation. Proceeding then as in Ref. p^ . 
one can verify that the divergent contribution to the velocity field, coming from the integration over large r, is also 
coordinate-independent, so that the condition c) is still met. Finally, it is not difficult to show that the expression 
([T^ cannot be further simplified following the lines of Ref. [l^l by omitting the additional potential contribution after 
the spatial differentiation: it turns out that this contribution satisfies the condition a) only in the steady case. We 
shall return to this point later in Sec. IV Al 



A. Definition of the vorticity mode 



Having obtained an explicit expression for the vortex component of the burnt gas velocity for a finite S, we have to 
consider the question of the transition to the limits A ^ oo, T oo. Generally, the rule these limits are to be taken 
depends on the problem under consideration. In the case of unsteadily propagating flames, this issue is complicated 
by the fact that the expression found for the vortex component is essentially nonlocal, both in space and time. The 
latter nonlocality shows itself explicitly through the r-integration in Eq. p^ . and is to be expected from the very 
outset. In fact, appearance of time non-locality is inevitable, because there must exist some mechanism transferring 
the influence of the flame history onto its current state. Such mechanism is unnecessary only in the stationary case, 
where all information about the flame past is in a sense left at the infinity downstream. Furthermore, we have seen 
that the time dependence of the function M(x,t) is partially transmuted into coordinate dependence, so the time 
non-locality entails naturally an essential spatial non-locality. This is again in contrast with the steady case, where it 
turned out to be possible to find a local on-shell expression for the vortex component [Cf . Eq. (37) in [l2| . 

Although the parameter T does not appear explicitly in Eq. p9|) . it should be kept in mind that in the course of 
derivation of this expression, the time dependence of the M(x,t) function has transmuted into dependence on the 
spatial coordinates. Hence, in order to preserve property c) of the potential component, the limit T oo is generally 
to be taken under assumption of vanishing of the function M{x,t) for t — oo. As to the limit A oo, any such 
assumption would be irrelevant because of the flame periodicity along the I'-axis. To ensure convergence in this case, 
we introduce an intermediate regularization of the x-integral replacing Eq. (fTO|) by 

+A T+ 

= W I die-^"-^ I dTM{i,t-T), (20) 



4-'^ 7 dxk 



-A 



where /i > is a sufficiently large parameter. We then take the limit A oo and define the vorticity mode as the 
analytic continuation of (j20p to the value /i = along the real axis in the complex /i-plane. Replacing also the symbol 
= by the equality sign, the deflnition thus reads 

+ 00 -r+ \ 

die'"''-^ J dTM{i,t-T)[ . (21) 




Let us show that the given deflnition using analytic continuation in fi respects the properties a) - c). Rewriting 
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Eq. ini) as 

+Ao r+ 



J dx-^ j dTM(x,t- t) 

dxk J 



— An +00 



-00 +Ao 




where ^0 > is arbitrary, and comparing with Eq. ()19p one sees that taking the hmit A ^ 00 followed by the analytic 
continuation in fi does not change the vorticity distribution in the arbitrarily large domain |a:| < ^o- Therefore, the 
above analytical operations amount to addition of some potential field, so that b) is met. Furthermore, this field is 
bounded in the sense of c). To see this, consider the quantity 



d f 
— / dTA/(i,t-r) 



(22) 



Using Eq. (jl5p and the Newton-Leibnitz formula, this expression is a combination of the functions Af (i, t~ t±) times 
spatial derivatives of t±. It follows from Eq. (fT51) that \t± | = r/i;+ — s- cx) for 00. Hence, if we assume that M (i, t) 

is exponentially bounded in a vicinity of the point i = 00 in the complex t-planc, i.e., \M{x,t)\ < e'^'*' for some c > 
and \t\ 00, then there exists a large enough /i such that the i-integral in Eq. (|2ip converges. On the other hand, 
the function M{x,t) is periodic with respect to x as the result of the flow periodicity. Therefore, all singularities of 
the expression in the curly brackets in Eq. (j2ip are off the real axis in the complex yu-plane, except possibly for a 
simple pole at ^ = 0. The latter corresponds to an additive constant, Bk, in the quantity (|22p . The appearance of 
such a term is not forbidden by the requirement of periodicity in x. Bk is independent of x by virtue of the same flow 
periodicity. But after the x-integration this term gives rise to a contribution of the form B^/iJ. + hk{x), where hk{x) 
vanishes for /i — !■ 0. Since Bk/ fJ. disappears upon spatial differentiation, Dv" can be continued analytically to /i = 0, 
so the property c) is fulfilled indeed. Finally, the divergence of the vorticity mode 

^ - J dTM{x,t- T) > 

LI = 0+ 




4 J r dxk J 



is proportional to /.i. On the other hand, one has r^/r = — sign(i)(5i.i + 0(l/|i|). Therefore, the only term contributed 
by the integral, that survives after continuation to /i = 0, is an a;-independent constant proportional to (SuSikBk/ i-i). 
Thus, div v" = const, and the condition a) is satisfied. 

IV. CLOSED DESCRIPTION OF NON-STATIONARY FLAMES 

We are now in position to write down an integro-differential equation relating the on-shell values of the fuel velocity 
and the front position function. Let us introduce the complex variable z = x+iy, and the complex velocity lu = u+iw. 
By virtue of the properties a), b) the complex quantity dio^/dz, where = + iw^ , is an analytical function of 
the complex variable z in the downstream region. In conjunction with the property c), analyticity of du^ /dz can be 
expressed in the form of the following dispersion relation [ll|, [l^] 

1 + (c^P)' = 0, (23) 

where the prime denotes ^-differentiation, and the action of the operator IK on an arbitrary function a{x) is defined 
by 



(via) {x) 



+ 00 

x-x^ i\][x,t) - J{x,t}\ 
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slash denoting the principle value of the integral. 3{ has properties similar to the Hilbert operator H (and is effectively 
the Hilbert transform along the front). In particular, it was proved by [T^ that 

3<^ = -l. (25) 

The identity (|23p relates in a complicated way the on-shell values of the burnt gas velocity and the flame front 
position. The fuel velocity also satisfies the conditions a) - c), this time in the upstream region. Indeed, a) is just 
the differentiated Eq. b) follows from the Thomson theorem and the boundary conditions upstream, and c) is 
true because the fuel velocity is bounded. The consequence of these properties is the following dispersion relation for 
oj- = U- + iw- 

(l -ijc) (cj_)' = 0. (26) 

Denote [v] the jump of the gas velocity across the flame front, [v] = v{x, f{x,t) + 0) — v{x, f{x,t) — 0). Then the 
sought equation for cu^ , / is obtained by substituting 



in Eq. (gSl), and using Eqs. (gl 



+ 00 Tj- 

2{uj^)' + (l + i3<) {[lu]-^ j diek-^ j dTM{S:,t-T)} =0, (27) 

— OO r_ 

where Cfc = £2fc+i£ifc) and we omit for brevity the regularizing factor e"'"' in the integrand as well as the accompanying 
symbol of analytic continuation. In the last term on the left, the argument y is understood to be set equal to ,f{x,t) 
after the spatial partial differentiation is performed, but before the ^-differentiation denoted by the prime. The value 
of vorticity at the front and the normal velocity of the burnt gas, entering the function M{x,t — r), as well as the 
velocity jumps at the front are known functionals of the on-shell fuel velocity [l^, [l^. For instance, for zero-thickness 
flame fronts one has 

vl - 0, , [w] = -/-^, (28) 



D d f f'\ d 



where 

_D _a / 

Dt^di^y^' ' N J dx 

Thus, the complex Eq. (|27p gives two equations for three fimctions w^{x,t), u_{x,t) and f{x,t). Together with the 
evolution equation 

(i3_ • n) = 1 + 5(m_ , w_ , /') , (30) 

where 5* is a known functional of its arguments, proportional to the flame front thickness (or. rather, the cut-off 
wavelength Ac), Eq. (j27p provides one with a closed description of unsteady flames in the most general form. Its 
application to various particular problems is given in the next section. 

Before proceeding, the following point is worth to comment on. In the analysis that led us to (|27p . the front was 
represented by a graph y = f{x, t), which excludes the overhangs or the fronts that double back on themselves. This 
over-restrictive assumption, adopted so far for simplicity, can be relaxed as follows. Let us parameterize the flame 
front by a real parameter ^ so that {x{^),y{^)) be a diffeomorphic mapping of the interval — oo < ^ < -|-oo onto the 
front at time instant t. Then, setting x{£,) + iy{^) ~ Z{(,,t), we define the new metric coefficient N{£^,t) in terms 
of the infinitesimal arclength along the front, dl, by dl = N{^,t)d£,; in other words, N{^,t) = \d^Z{^,t)\ , which is 
nowhere singular for smooth fronts. As long as the fuel and the burnt-gas regions remain connected domains, the 
Cauchy theorem guarantees the existence of a generalized operator M such that (1 ± i'K)duj±/d^ ~ 0. Specifically, 
when acting on a smooth a{^,t), the new !K produces 



(31) 



10 



instead of Accordingly, if M(^, t) is still defined as M = N{^, t) cr+{C, t) (C, t), Eq. ([271) is formally unchanged, 
except that the prime now denotes d/d^. Of course, the time derivatives, now at fixed ^, must be handled in a way 
consistent with the new representation of the flame front. Yet such a re-parameterization does not capture situations 
when isolated pockets of unburnt fuel form, because the fresh domain then ceases to be path-connected. Unfortunately, 
as long as a local propagation law [Cf. (|30p] is employed such a phenomenon cannot be excluded a priori: in no way 
can a flame element "know" that another one is to produce a "head-on" collision. 

The above re-parameterization is unnecessary for the wrinkled fronts considered below. 



V. EQUATIONS (HH), jSOD IN LIMITING CASES 



To give a consistency check for Eq. ([27|) and also to gain a deeper insight into the structure of this equation, we 
use it below to derive anew classical results on flame front dynamics. 



A. Darrieus-Landau instability of zero-thickness flames 

Let us consider first the classical linear stability problem of zero-thickness planar flame propagation 0, [1] . In this 
case, = 9, N = 1, while the linearized on-shell vorticity ([29|) reads 



Accordingly, expression (|2ip simplifies to 



9~ \ ~dr ^ dtdx ) 



Since the integrand here is a first order quantity, it is sufficient to calculate t± for a plane front which is assumed to 
be at 2/ = 

y , . \x-x\ 

r± = ^±.^. (33) 

Let the disturbance be periodic in x, growing exponentially with time. Spatial periodicity of the linear problem is 
most conveniently represented in the complex form, in which case 

fix, t),u^ (x, t), (x, t) ^ 6^'=-+''* . (34) 

However, it should be kept in mind that the coefficients in Eq. (|27|) are also complex. In order to preserve the right 
complex structure of this equation, all the functions involved are to be written in the form containing no imaginary 
coefficients. To find w^, we have to evaluate the following integrals 

+OC y/9+i\x~x\/9 

d 



y/ 9—i\x~x\/ 9 

for > and A; — 1, 2. A straightforward calculation gives 



y\I^L + ik + lu I a —i^L + ik + iv/tiJ 

y\I^L + ik + IV I u —^L + ik + iv/OJ J 

where {v — s- ~u) is shorthand for the preceding parenthesis with v changed into its opposite. The on-shell values of 
these functions analytically continued to /i = are 

/r(x,0,M) = /2(x,Q,t,0) = ^e™^^^{§^. (35) 
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Using this in Eq. (132]) yields 

ik (e^i) 



To put these expressions in an explicitly real form, it is sufficient to write 



[vw- + ikvf) , (36) 
{vw^ + ifci//) . (37) 



Next, we note that in the linear approximation, the operator becomes just the Hilbert operator whose action 
on the harmonic functions is defined by 

H exp{ikx) — ix(k) cxp(ifca;) , k ^ , (38) 

where 

x(fc) = 

Equation ((77|) thus becomes 




2(._)-+(l + .i^)^M + V 27;2^ -0. (39) 



Extracting the real and imaginary parts of this equation, we find 



(w_)' = H(u„)'. (41) 
Finally, upon linearization the jump conditions (j28p simplify to 

M = - 1 , N = -/'(^ - 1) , (42) 
while the linearized evolution equation reads 

..-|.l. (43) 

Inserting these into Eq. (|40|) and then substituting (7+ = ^(^^ ^) / (^{'ix{k)i''^ f + iki' f) leads after some simple algebra 
to the equation 

t^„^ + 2i^\k\-{e-l)k^ = Q, (44) 



which is nothing but the famous Darrieus-Landau dispersion relation determining the perturbation growth rate as a 
function of the wave number and the gas expansion coefficient d, 0] 



f^?-i-lj|fc|. (45) 

The effects related to the finite front thickness can be explicitly accounted for in the above linear analysis by including 
terms linear in /" in the right hand side of Eqs. This would changed the last term in Eq. by the extra 

factor (1 — |fc| Ac/27r), and make v roughly parabolic for all |fc| < 27r/Ac; max(i^) enables one to estimate the typical 
growth time as t^L — 2\c/n{V0 — 1), to be used later (Sec. 6). 
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Next, let us return to the remark made after Eq. p4|) . It was mentioned there that even if one is interested only 
in the on-shell values of the vorticity component, the r-dependence of the function M(x, t — r) cannot be neglected. 
Evidently, doing so amounts simply to rewriting the equation obtained by Kazakov [lll.[l3| for steady flames in terms 
of the local current flow velocity relative to the front. It is not difficult to verify that in the case under consideration, 
this would change the coefficient of i'^ in Eq. (|44)) to the wrong value [6 —l)/9 . This change is thus a reflection of the 
memory effects encoded in the function M{x, t — t). Equation ([27|) properly takes into account these effects, correctly 
reproducing the Darricus-Landau relation, and automatically captures all the aspects of flame dynamics related to 
inertia. 

At this stage of the analysis it is appropriate to pause and discuss the meaning of the analytic continuation appearing 
in the definition (|2ip in somewhat more detail. This continuation is used to make the i-integral meaningful in the 
limit A —t OO. One can avoid using this analytic means if the formal result of improper integration along the infinite 
flame front is treated in the sense of distributions. Indeed, in this case expressions (j35|) for the vortex component 
would contain additional terms proportional to 5{k + v /9) or 5{k — v jff) coming from the integration over large is. 
To see that these terms are inconsequential we recall that the above consideration of the single fc-mode evolution 
is not completely adequate from the physical point of view, because in practice one always deals with wave packets 
consisting of a continuum of wave numbers fc, rather than with a single mode. This means that the physical expression 
for the flame front position with the given fcg is obtained by integrating the found solution f{x,t) with some weight 
over a small but finite range Ak around fco. Upon this integration all terms involving d{k ± v/O) disappear, because 
V = ±k9 are not roots of the Darrieus-Landau relation. The two procedures are thus equivalent. 



B. The small {9 — 1) expansion 



Let us next consider the case of small gas expansion. We will verify that within the framework of the asymptotic 
expansion with res2ect to 6' — 1 = a, Eq. (|27|) reduces at the first post-Sivashinsky order to the well-known Sivashinsky- 
Clavin equation [ij]^. To perform the asymptotic expansion we recall that the cutoff wavelength for the short 
wavelength perturbations Ac is of the order l/a. This means that the wave numbers involved are 0{a). In other words, 
spatial differentiation of a flow variable raises its order by one; in particular, /' = 0{a). Also, Eq. (|45p tells us that for 
small a, the perturbation growth rate v = ka/2 = 0{a^), so the order of a flow variable is raised by two upon time 
differentiation. It follows then from Eq. ([30|) [with S* = 0] that u- = 1 + 0{a^), while potentiality of the upstream 
flow implies that w- = 0{a^) [this is clearly seen from the dispersion relation (dH)]. The flrst post-Sivashinsky 
approximation corresponds to retaining terms of the fourth order in Eq. (j27p . or equivalently, 0(a'^)-terms before the 
spatial differentiation. It was shown by [l3| that to this order, 5{ becomes just the Hilbert operator. Therefore, like 
in the linear case considered above, the real and imaginary parts of Eq. (|27p arc readily separated to give 

[u] - H[w] + {eikH - £2fc)^ J di-^ J dTM{i, t t) \ = , (46) 

— OCT- ) 

= H{u_)' . (47) 

Since the quantity Jk = dk J^^ drM is under the sign of spatial integral, we need the function M[x, t) to the fourth 
order in a. To this accuracy, it is equal to the on-shell vorticity 

M = o+ = -{6-l) [f' + f'f") , 



^ It is worth noting that the small expansion parameter used by Sivashinsky and Clavin is 7 = {9 — 1) /9 = a/ {1 + a), rather than a. The 
reason for switching from a to 7, found sometimes in the literature, is that the latter would improve the expansion accuracy/convergence. 
The argument given in this connection, namely, that 7 < 1 for all 9, while a is small only for 9 close to 1, is not quite correct. Validity 
of an expansion is determined not by the value of the expansion parameter, but by the relative value of the terms neglected in the 
course of the calculation. It is true that expanding the decreasing function f{9) = 1/9 = 1/(1 + a) in terms of «/(l -I- a) instead of a 
is an improvement. However, it is not in the case of the increasing function f[9) = 9 — 1 = a, which certainly appears in the governing 
equations [see Eqs. H27| l. I|28| l] . Thus, whether a change of the expansion parameter improves the expansion accuracy /convergence is a 
question of the structure of the whole perturbation series, which cannot be resolved from the knowledge of its first few terms. It is one 
of the goals of our approach to make questions of this kind accessible for theoretical analysis. 
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the dot denoting differentiation with respect to the time t. As this quantity is already 0{a'^), t± can be talcen in the 
form (|33p with 9 — 1. Thus, performing the spatial differentiation and setting y = gives 

Ji = ixix — x) {a+{x,t — i\x — x\) + a+{x,t + i\x — x\)} , 
J2 = {(T+{x,t — i\x — x\) — (j+{x,t + i\x — x\)} . 

Let us show that imaginary parts in the argument of a+ can be omitted. Note that within the asymptotic expansion 
in a, dependence of cr+ on i\x — x\ can be treated perturbatively. Indeed, let us write formally 

cr+{x, t — i\x — x\) = (J+{x, t) — i\x — x\&+{x, t) — —{x — i)^(7+(i, t) + ■ ■ ■ . 

To assess the relative order of consecutive terms in this series, we have first to get rid of the factors containing explicit 
coordinate dependence, which can be done by successive integration by parts with respect to x. One has, for example, 

+00 +00 fx \ 

dx\x — x\a j^ix = d, \ dxia+{xi,t) \ \x — x\, 



+00 / i ^ 

J dxx{x — x) i J dxia+{xi,t) 



where according to the discussion in the preceding section, the contributions of infinitely remote parts of the front 
have been omitted. Since each time differentiation adds two powers of a, while spatial integration subtracts only one, 
we conclude that the above expansion for (7+ is effectively an asymptotic series in powers of a. Hence, to the fourth 
order in a 



Ji — 2ix{x — x)(T+{x,t) , J2 = , 



and so 



+ QO -\-00 

diJi = -2i{9 - 1 



I dixix^i) {f' + f'f") {x,t) 



-00 
+ CXD 



= -4^(^? - 1) y dxSix -£)(^f + ^ j (i, t) 



= -i^{0-l) (/ + ^) (^'*): 

where the boundary terms at i = ±00 were omitted again in the course of the integration by parts, and the relation 
dx{x — x)/dx = —25{x — x) was used. Inserting this together with the velocity jumps into Eqs. ((^ . (gT]) yields 
within the given accuracy 



2u- + {e-l)[Hf' + f)=Ci{t), (48) 
2w- + {e-l)(-f' + Hf)=C2{t), (49) 



where Ci_2(0 are two "integration constants" (some x-independent functions of time). Since the left hand side of ((49 
is odd in x, one has C2 = 0. Using these formulas in the evolution equation (|30p written in the form 



f'2 

- / - f'w- = 1 + V 



leads to the equation for the function f{x, t) 



'-i^f+'^r + Cit) = -'-^Hf', (50) 
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where C{t) = 1 — Ci(i)/2. The function C{t) can be found by averaging the obtained equation along the flame front: 

Equation (jSO]) is nothing but the Sivashinsky-Clavin equation [l3|, with the term C{t) added according to [20| . 
Account of the transport processes inside the front would have added a term proportional to Ac/" to the right hand 
side of Eq. ([50l) . 

Finally, it is interesting to note that the memory effects are insignificant not only at the lowest order of the small 
{9 — l)-expansion, but also at the first post-Sivashinsky approximation considered here. Namely, a direct calculation 
shows that replacing M{x,t — r) by M{x,t) in Eq. (|27|) does not change Eq. ([50]) . This is natural, because the 
condition a <C 1 implies slow dynamics. However, this is not the case already in the next order, which is possibly one 
reason why the memory effects are often overlooked. 

C. Flame propagation in time-dependent gravitational field 

Depending on the direction of flame propagation, gravitational field leads either to strengthening or damping of the 
Darrieus-Landau instability. Also, the influence of sound waves on the front dynamics can be described effectively as 
the flame propagation in a time-dependent gravitational field [Tsj . Let g{t) denote its strength, with the convention 
that g > corresponds to a stabilizing effect [see Eq. ^S^]. Inclusion of the gravitational field does not change the 
flow equations so that their consequence, Eq. (P7)) . has the same structure. The jump conditions for the 

gas velocity are also left intact by gravity. The only place where g{t) appears in our approach is the expression of the 
on-shell vorticity (as a result of baroclinic effects inside the front). Namely, the gravity-induced jump in this quantity 
to be added to the right hand side of Eq. ^ is [HI'* 

Aa+ = -t^.9W/'(x,t). (51) 

If development of the Darrieus-Landau instability is suppressed by the gravitational field, then it is natural to consider 
harmonic front perturbations instead of the exponentially growing ones used in Sec. IV Al The choice between the 
two representations is just a matter of convenience. Namely, it will be seen below that although the intermediate 
procedure of analytic continuation in fi is somewhat different for v imaginary, the final equations for the front position 
are simply analytic continuations of one another with respect to frequency. 

Because of the time dependence of g{t), the function f{x,t) is now to be taken as a superposition of an arbitrary 
number of harmonics 

-|-oo 

ikx — iut 



fix,t)= j dujdkf{k,uj)e"'''-''^\ (52) 

— OO 

It is also convenient to introduce the Fourier decomposition of the function G{x^t) ~ g{t)f'{x^t) 

G{x,t)^ I dujdkG{k,uj)e''''=~"^K (53) 



To perform the analytic continuation with respect to ^ in Eq. (|2ip . we will assume that f{k,iu), G{k,Lu), considered 
as functions of w, vanish outside a large but finite frequency band ^ ujq. Then, choosing /j, > uJo/9, inserting 
the above Fourier decompositions in the linearized expression for the on-shell vorticity, performing i-integration in 
Eq. (|2T|) as before, and continuing the result analytically to ^ = 0, we find 



W; = —I ■ 



^ ^ f dLudke^'^''-^-'-^^^^^{i-ixik)uj' + kLu)fik,Lu) + Gik,Lu)} , 



-OO 



+ 00 

I dLodk e"— * ^ ^^^^^ { i-ix{kW + kLo)f{k, Lo) + G(fc, u) } 



More generally, Acr^ = ^^njj^ {r(x, /{x, t), t)}' , where r[x,y,t) is the gravitational potential, g = — VF. 
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The first terms on the right of these equations are just analytic continuations of the expressions ((36)) . ([37]) to the 
imaginary value of the growth rate: v ^iuj. As in the Darrieus-Landau problem, Eq. (j27p splits into two real 
equations 

2(w_)' + {M - H[w] -ul + Hwl^' , = H{u_y , 

which after substitution of the found expressions for together with Eqs. (H^ . yield 

1 



+ 00 



dtodk e 



ikx — iujt 



-u? + 2iuj\k\ + {e - i)k'^ 

+ ^—r^ikx{k)Gik,Lu] 



kf{k,uj) 



ik 



0. 



Acting on this equation by the operator 



id 9 , , , 



and taking into account definitions (|52p . (j53p . we arrive at the well-known equation for the front position [IE 



-f-2Hf' + {6-l)f" 



-g{t)Hf' = C{t) , 



(54) 

where C(i) is a function of time appearing as the result of the spatial integration symbolized by ^ dx. As before, 
C{t) can be found by averaging Eq. ((54)) along the front. Equation (l54l) is well-known to be the key ingredient in 
understanding the parametric flame response to oscillating g{t)s jiol Il6l|. in which case the inertia (hence, memory) 
effects play the central role. 



D. Steady flame propagation 

In the case of stationary flames, the A/-function is time-independent, whence the r-integration in Eq. (|27[) is trivial. 
One has 

-I-CX3 ''"+ -(-OO 



where 



I dxcfc— dTM{x,t~T) = -- / dx— ^efc — Jr2 - ^ ^ 

dxk J 2 J v+ dxkV < 

r_ — oo 

-t-oo 

^ dxM{x) ^'^^ , 



,3, ^ [in n^J±\ . (55) 



This vector satisfies 



i.e., (3 is the unit vector orthogonal to v+ . In addition to that, f3 changes its sign at the point x satisfying Eq. ([T 
Hence, its scalar product with the complex vector can be written on-shell as 

ekPk = -—x{x - x) . 

Differentiation with respect to x using the formula x'ix) = 2S{x) makes the i-integration trivial, so that Eq. (j27p 
takes the form 

2(.._)'+(i + .^c)|h'-^^^^|=0 (56) 
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which is exactly Eq. (42) of Ref. [l3|. 

In connection with Eq. (|56p one can indicate why viscous effects have virtually no influence on the shape and the 
speed of steady flames, except possibly a small one on the cut-off wavelength (Ac) itself, the only internal reference 
length of our problem. The Reynolds number based on Ac and the burnt gas properties and speeds is normally about 
10^. The viscous length Avis - the distance at which the vorticity present behind the front-crests of transverse size 
~ Ac dissipates noticeably - thus exceeds Ac by two orders. This is definitely too late (in terms of the Lagrangian 
time r) to modify Eq. (|56p . whose validity only requires that cr (T-(_ in an infinitesimal layer (r — > 0+) downstream 
of the front (r = 0). As to unsteady flames the picture is more involved, even though the potential fuel flow is still 
unaffected directly by viscous effects. The problem would certainly deserve further studies, for one cannot a priori 
exclude viscosity- related large-scale phenomena occurring far from the front, that could nevertheless influence its 
dynamics indirectly. We merely mention here that the spontaneous flame dynamics is initially little affected, because 
the Darrieus-Landau time tuL (see Sec. 5.1) is well shorter than ivis = \is/0: the decay of M{x,t — r) for r ::g> 1 is 
fully controlled by the wrinkle growth itself. 



VI. STABILITY ANALYSIS OF GENERAL STEADY FLAME PATTERNS 

Let us use Eq. (|27p to derive an equation governing propagation of small disturbances along a given steady flame 
pattern. As in the Darrieus-Landau or stationary problems considered above, Eq. (|27p essentially simplifies in this 
case, because the time non-locality is no longer a complication. Indeed, since there is no external time- varying field (as 
the stationary regime is assumed to exist) and the flame disturbance is small, it is sufficient to consider perturbations 
having the form 

Sfix,t)^f{x)e''\ Sw-{x,t) = w{x)e''\ 6u^{x,t) =: u{x)e''\ 



in which case the time dependence is prescribed, and the r-integration in Eq. (I27|) is readily done. Let us gather the 
functions f{x,t), W-{x,t), U-{x,t) into an array {£^a{x,t)}, a = 1,2,3: 

(6,6,6) = (/,w_,u_) , 

and denote ^i^-* (x) the given solution of Eq. ([SS]) . We have to linearize the left hand side of Eq. (P7|) with respect to 
S^{x,t) = £^{x)e'^* = {f ,w,u)e'^* . Depending on whether the function M{x,t) is varied, one finds contributions of two 
types. If 6 comes from terms other than M, then M(i, t) is evaluated for <^ = S}-'^'^ and hence is time-independent. As 
we saw in the preceding section, Eq. (|27p simplifies in this case to Eq. (j56p . so the corresponding contribution to the 
variation is given by the left hand side of Eq. (|56p linearized around ^^"^ , with M kept fixed. The other contribution 
comes from variation of M{x, t) and has the form 



where Maix) is the differential operator obtained by linearizing the function M{x,t) around the stationary solution, 
and changing d/dt —f v afterwards. A straightforward calculation gives 



+ 00 



Am = (l + i?i:j y diM„(£)Cc«(£)^ |exp e-*^j x(-^-i) 

— oo 

where e tt, +7r] is the angle between the vectors r, defined positive if the rotation from Vj^ to r is clockwise. 
It is not difficult to verify that the x-differentiation of the step function in the latter expression gives rise to the term 
which is just the variation of the left hand side of Eq. (|56p under variation of the function M{x, t). Furthermore, one 
has 

. d{rsm<j)) 
(rsm0) = Ti = TiP.xix - X) , 

OXi 

where = Sikrik is the unit vector tangential to the fiame front. Taking into account also that 

SikPi = X(X - X) , 

V+ 
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wc find 

(rsin0) = eikUkPzXix - x) = = — 

v+ v+ 

Similarly, 

, Vk+Tk _ V+ 



(r COS (/))' 



v+ v+ 



Putting all these results together and omitting the factor e''*, we obtain the following equation for the x-dependent 
parts of the perturbations 



+00 



= ^fl + z^c) / diM^{i)Ui) + 7+^"^+ exp ( --e-'^ ) x(x - ~x) , (57) 



where (5/i5Cq(-^^) denotes the functional differentiation, and a) = u^iw . Together with the linearized evolution equation, 
Eq. (j57p can be used to carry out stability analysis of general steady flame configurations, given by the solutions of 
Eqs. (|5D|) . (|55|) . To the best of our knowledge, no such closed equation as ([57)1 has yet been derived to handle the 
problem. Its applications will be presented elsewhere. 



VII. DISCUSSION AND CONCLUSIONS 



The results presented in this paper solve the problem of non-perturbative description of unsteady premixed flame 
propagation with arbitrary gas expansion. Supplemented by the evolution equation, Eq. (P7|) gives the closed de- 
scription of unsteady flames in the most general form. Thus, as in the stationary 2D case, the dilemma mentioned 
in Introduction is resolved in negative for 2D unsteady flame propagation. The conclusion that the detailed bulk 
structure of the gas flow is actually unnecessary for describing front dynamics is even more striking in the latter 
case. Indeed, the highly complicated vortex flow downstream continuously changing in time is naturally expected 
to have an exceedingly complicated nonlocal influence on the flame front evolution. We proved, however, that all 
necessary information about this influence is encoded in the complex history of the combination M = Nv^<j+ . In 
this connection, a curious circumstance is worth mentioning. As we saw in Sec. IIIIl the vortex component depends on 
spatial coordinates through the complex combinations t± appearing as the limits of integration in the complex time 
plane. It is easy to see that for any curved flame configuration, there are always points at the front, where < 0. 
For such points, the time argument of the function M{x, t^r) in the integrand of Eq. (j2ip has real part > t. In other 
words, integration over such points is in a sense looking into the flame future. This does not lead to any conflict 
with causality, because the corresponding contribution is eventually annihilated by the operator (1 + i'H) in Eq. (j27p . 
Indeed, the vortical part of v" comes from integration over x such that the vectors r and ■u^. are almost parallel, i.e. 

w 1, and therefore, the contribution of points with < is a pure potential satisfying Eq. (|23p . Retaining this 
contribution in intermediate formulas is necessary to guarantee continuity of the potential component. 

The paradox as to the seemingly "teleological" structure of (^7)) is closely related to the analyticity properties of 
M{x, t) in the complex i-plane, that allowed us to simplify Eq. (|14p to Eq. (US]). Our considerations were carried under 
the very weak assumption of exponential boundedness of this function, which is certainly sufficient for investigation of 
any flame propagation phenomena. In particular, the linear stability problem analyzed in Sec. IV Al gives an example 
of M{x,t) which is analytic in the complex plane, so that the question whether the contour deformation is legitimate 
does not arise at all. In this simplest case knowing the exponentials exp(j/t) in Eq. (|34p at current time t allows one 
to predict their future, so that the "teleological" question does not arise either. Things are more interesting, however, 
in the case of flame propagation in time-dependent gravitational field. Suppose that the experimentalist plans to 
leave the burner fall freely at some time instant tQ. This means that g{t), and hence, M{x,t) [Cf. Eqs. pH]) . ([ST]) ] will 
have singularities at tq = to ± * At, where At is of the order of duration of the field switching off. As is shown in the 
Appendix, crossing these singularities by the contour C_ U C+ in Fig. 2 gives rise to a potential contribution that 
satisfies the conditions a)-c) of Sec. Illli and hence does not change Eq. (|?7|) . thus resolving the causality issue. 

To conclude, Eq. (j27p opens a wealth of key developments in theoretical and numerical combustion (extended 
propagation laws, coupling with acoustics, burners, etc.), not to mention the other fronts evoked at the beginning 
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of this paper. In particular, Eq. (j57p allows direct analytical investigation of small disturbances propagating on 
steady front patterns such as Bunsen- or V-flames in 2-D configurations. Extension of the above results to the three- 
dimensional problems is still an open question, one of the main difficulties being the generalization of the (1 + iOi) 
operator projecting out the potential contributions of the burnt-gas flow. 
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APPENDIX A 



When deriving the expression (j2ip for the vortex mode in Sec. IIIII we have omitted the contribution of the contour 
integral in Eq. (|17p , retaining only the singular contribution of the poles t± . That this operation respects the property 
c) was already shown in Sec. IIII Al We will now prove that Eq. ([21]) does reproduce correctly the near-the-front 
distribution of vorticity of the burnt gas flow. In essence, the subsequent calculation reproduces the consistency check 
given in the appendix A of • First of all, taking into account the formula 



dxj 



V+ \ V+ 

which is readily verified using the definitions (fT5|) . ([55]) . one has 



(Al) 



SkidkVi = -T— 




M{i,t-T+) { A-z^ 
v+ 



+M{i,t-T^) [(3, + i^ 



.Vi+ 



■1=0+ 



f+CX3 
f e"^"" / V 

/ di__M(£,t-r+) A-i- 
J v+ \ V 



Following the argument given in Sec. IIIII one performs the differentiation with respect to Xi under the sign of the 
i-integral, and sees that the derivative of e~'"" leads to an integral proportional to ^. Therefore, the only non-zero 
term contributed by this integral after ^ is continued to zero is an inessential a;-independent constant that falls off 
from the expression miDv'" . Next, it is easily checked that 

A--*^) =0, (A2) 
so differentiation of M{x,t— t+) gives zero, too. Note also that f3 can be written as 

A = ^^^i^x(eimnwm+) 7 (A3) 
since it is orthogonal to and changes sign at the point satisfying Eq. (|12p . Thus, we find 

SkidkVi = RB'I / dx^^M{x,t- T+)eikVk+£inVn+S{ehnriV.m+) 



vl 



+ 00 

dxM{x, t - r/v+)S{eimriVm+) ■ (A4) 
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The factor e"'"' and the symbol of analytic continuation have been omitted in the last expression because it is explicitly 
finite. The argument of the (5-function turns into zero when the vectors and are parallel. Near this point, one 
has approximately 

ehnriVm+ = rv+cj) . 

On the other hand, a simple geometric consideration shows that dx near the same point can be written as (see 
Fig. El) 

rv+dd) 

dx -- 



where all quantities are taken at the time instant t. Substituting these expressions into Eq. (|A4p . and taking into 
account the relation 



6{ax) — - — tS(x) 
a 



yields 



[ d^Si^)Mpl^im . (A5) 



If the observation point is taken at the flame front, r turns into zero together with (j). Recalling also the definition 
of the function M{x,t), we arrive finally at the desired identity 

{ekidkv^)^ = (7+{x,t) . 

Let us now return to the question of possibility to perform the contour deformation in the complex r-plane, used 
in the derivation of Eq. p^ . We note, first of all, that this question is concerned entirely with the structure of the 
potential component of the gas velocity. Indeed, under our general assumption of existence of a short wavelength 
cutoff, all functions involved are smooth functions of time (because any structure of finite size takes finite time to 
develop). Hence, these functions (in particular, the function M(i,t)) are analytic in a vicinity of the real axis in the 
complex i-plane. On the other hand, we know that the value of vorticity in any given point [x, y) near the front is 
equal to its value at the point on the front, satisfying f2 = 1, in which case the poles t± take the real value r /v^. For x 
in a vicinity of that point, t± belong to the analyticity domain of M (i, t), and hence, of the function !M(5;, r, t). Thus, 
the use of the Cauchy theorem and the contour deformation performed in Sec. IIIII are legitimate for these i, while 
integration over all other x gives rise to a potential contribution. This, however, does not conclude consideration, 
because the argument just given proves potentiality of this contribution in some vicinity of the given observation 
point. To prove it for all x £ [0, 1], the contour C_ U C+ in Fig. [2] should be moved to the left of t± for all x. Let 
us show that this is possible indeed under the assumption used already in the derivation of Eq. (j2ip . namely, that 
M{x,t) considered as the function of the complex t is exponentially bounded near t ~ oo. Suppose, for instance, 
that this function is meromorphic, i.e., has only poles of arbitrary order in the f-plane. Any pole of the order n in 
the function M{x,t) becomes an {n — l)th order pole with respect to r in M(i,T, t). Hence, crossing these poles by 
C- U C+ does not change the right hand side of Eq. p7)) unless n = 2 or n = 1. Consider first the case n = 2. Then 
M(i, r, t) has simple poles at some point tq and its complex conjugate Tq (tq may depend on i, t, but we do not write 
this dependence explicitly, for brevity). We are to show that their contribution to the right hand side of Eq. (|17p . 
given by the Cauchy theorem as ("c.c." stands for complex conjugate) 

— resM(a;,ro,tK \ ^ + c.c. , A6 

47r 1^ To - r+ To - r_ J 

gives rise to a field Vi that satisfies conditions a) - c). Substituting this expression into Eq. pT|) yields (we do not 
introduce intermediate regularization because the i-integral will be shown to converge) 

+A 

V^^~^-£ik / resM(£,To,t)afe i \ I + c.c. (A7) 

4 i L To - T+ To - T_ J 

-A 
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The property a) is evidently satisfied. To prove b) we write, using Eqs. (lAip . (|A2p and (jA3|) 

+A 

f r 1 1 

dk / resM(a;, To, t)9fe 



-A 

+A 



dx resM(i;, tq, t) 



To - T+ To - T_ 



-A 

+A 



(to - T+)2 (to - T_)2 
1 



2i / diresM(i,To,t)(5(ei,„ntim+) 



(to-T+)2 (to-T_)2 



since t+ = t_ when the argument of the dcha function is zero. Thus, rotV = . Last, M(x, t, t) is periodic in i, 
and therefore, so is its pole. Taking into account also that t± = 0{\x\), diT± = 0(1) for oo, one sees that the 

i-integral in Eq. (|A7[) is convergent in the limit A ^ oo for all x. 

In the case n = 1 the function M(i,T, contains a logarithmic singularity of the form a(i,t)ln(T — tq) (and its 
complex conjugate). Crossing this singularity leads to the 2tt jump in arg(ln(-)) for all points of the contour C_ U C+, 
located at one side of the point to. Hence, expression (jA6[) is replaced in this case by the following 



To 

2-!Tia{x,t) /" , f 1 1 



4tt 



dr i \ } + c.c. 



The proof of the properties a)-c) is exactly the same as before. 

Let us finally consider the case when the function M{x,t) has branch singularities. If these singularities are 
connected by a number of cuts so that M(x,t) is meromorphic in the cut T-plane, then so is the function M(i:,T, f), 
and moving the contour of integration beyond a cut results in a contribution to the right hand side of Eq. p7)) of the 
form 



1 

47r 



dT[M]{x,T,t) 



Co 



where [M](i,T, i) denotes the jump^ of the function M(5;,t, i) across the cut Cq. If this cut has finite length, then 
the above considerations again apply literally. However, in the case of an infinite cut, the T-integral is apparently 
divergent. This means that such cuts, if any, are only allowed in regions where M(i,T, t) satisfies more restrictive 
conditions than the exponential boundedness. We do not pursue details here, because physical significance of such 
cuts is not clear. 

Thus, the function M{x, t) is allowed to have any number of branch singularities in the complex t-plane, connected 
by cuts of finite length, as well as any number of poles of arbitrary order to justify the contour deformation used in 

secnni 
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FIG. 2: Extraction of 



the singularity in Eq. 



(|14p by contour deformation in the complex r-plane. 




FIG. 3; Near-the-front structure of the flow downstream. 



